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Abstract 

We develop a qualitative geometric approach to swimming at low 
Reynolds number which avoids solving differential equations and uses 
instead landscape figures of two notions of curvatures: The swimming 
curvature and the curvature derived from dissipation. This approach 
gives complete information for swimmers that swim on a line without 
rotations and gives the main qualitative features for general swimmers 
that can also rotate. We illustrate this approach for a symmetric 
version of Purcell's swimmer which we solve by elementary analytical 
means within slender body theory. We then apply the theory to derive 
the basic qualitative properties of Purcell's swimmer. 

1 Introduction 

Micro-swimmers are of general interest lately, motivated by both engineering 
and biological problems [9j El [161 El El El HI [HI HZ] • They can be remark- 
ably subtle as was illustrated by E. M. Purcell in his famous talk on "Life at 
low Reynolds numbers" [13] where he introduced a deceptively simple swim- 
mer shown in Fig. [TJ Purcell asked "What will determine the direction this 
swimmer will swim?" This simple looking question took 15 years to answer: 
Koehler, Becker and Stone [3] found that the direction of swimming depends, 
among other things, on the stroke's amplitude: Increasing the amplitudes of 



certain small strokes that propagate the swimmer to the right result in prop- 
agation to the left. This shows that even simple qualitative aspects of low 
Reynolds number swimming can be quite un-intuitive. 

Purcell's swimmer made of three slender rods can be readily analyzed nu- 
merically by solving three coupled, non-linear, first order, differential equa- 
tions [TB]. However, at present there appears to be no general method that 
can be used to gain direct qualitative insight into the properties of the solu- 
tions of these equations. 

Our first aim here is to to describe a geometric approach which allows one 
to describe the qualitative features of the solution of the swimming differen- 
tial equations without actually solving them. Our tools are geometric. The 
first tool is the notion of curvature borrowed from non-Abelian gauge theory 
|18j . This curvature can be represented graphically by landscape diagrams 
such as Figs. [3|5p 1 which capture the qualitative properties of general swim- 
ming strokes. We have taken care not to assume any pre-existing knowledge 
about gauge theory on part of the reader. Rather, we have attempted to use 
swimming as a natural setting where one can build and develop a picture of 
the notions of non-Abelian gauge fields. Purcell's original question, "What 
will determine the direction this swimmer will swim?" can often be answered 
by simply looking at such landscape pictures. 

Our second tool is a notion of metric and curvature associated with the 
dissipation. The "dissipation curvature" can be described as a landscape di- 
agram and it gives information on the geometry of "shape space" . This gives 
us useful geometric tools that give qualitative information on the solutions 
of rather complicated differential equations. 

We begin by illustrating these geometric methods for the symmetric ver- 
sion of Purcell's swimmer, shown in Fig. [2j Symmetry protects the swimmer 
against rotations so it can only swim on a straight line. This makes it sim- 
ple to analyze by elementary analytical means. In particular, it is possible 
to predict, using the landscape portraits of the swimming curvature Fig. |3j 
which way it will swim. In this (Abelian) case the swimming curvature gives 
full quantitative information on the swimmer. We then turn to the non- 
Abelian case of the usual Purcell's swimmer which can also rotate. There 
are now several notions of swimming curvatures: The rotation curvature and 
the two translation curvature. The translation curvatures are non-Abelian. 
This means that they give precise information of small strokes but this in- 
formation can not be integrated to learn about large strokes. This can be 
viewed as a failure of Stokes integration theorem. Nevertheless, as we shall 
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explain, they do give lots of qualitative information about large strokes as 
well. 




Figure 1: Purcell's swimmer. The swimmer controls the angles 16^2 1 < 7T. 
The location of the swimmer in the plane is determined by two position 
coordinates x 1,2 and one orientation x° = <fi. The lengths of each arm is i 
and the length of the body £q. 



2 The symmetric Purcell's swimmer 

Purcell's swimmer, which was invented as "The simplest animal that can 
swim that way" |13| . is not simple to analyze. A variant of it that is simple 
to analyze is shown in the Fig [2j The swimmer has four arms, each of length 
£ > and one body arm of length Eq, (possibly of zero length). The swimmer 
can control the angles 6j and the arms are not allowed to touch. Both angles 
increase in the counterclockwise direction, < 6j < tt. 

Being symmetric, this swimmer can not rotate and can swim only in 
the "body" direction. It falls into the class of "simple swimmers" which 
includes the "three linked spheres" of Najafi and Glolestanian [9] and the 
Pushmepullyou p], whose hydrodynamics is elementary because they can 
not turn. 

Let us first address Purcell's question "What will determine the direction 
this swimmer will swim?" for the stroke shown in Fig. [3j In the stroke, 
the swimmer moves the two arms backwards together and then bring them 
forward one by one. The first half of the cycle pushes the swimmer forward 
and the second half pulls it back. Which half of the cycle wins? 
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Figure 2: The symmetrized Purcell's swimmer can not rotate and can only 
move along the horizontal body-axis. The length of each arm is i and the 
length of the body £q. The position of the swimmer is denoted by X. 

To answer that, one needs to remember that swimming at low Reynolds 
numbers relies more on effective anchors than on good propellers. Since one 
needs twice the force to drag a rod transversally than to drag it along its axis 
|10j . an open arm 9j ~ ir/2 acts like an anchor. This has the consequence that 
rowing with both arms, in the same direction and in phase, is less effective 
than bringing them back out of phase. The stroke actually swims backwards. 
This reasoning also shows that the swimmer is Sisyphian: it performs a lot 
of forward and backward motion for little net gair^] 

The swimming equation 

The swimming equation at low Reynolds number is the requirement that the 
total force (and torque) on the swimmer is zero. The total force (and torque) 
is the sum of the forces (and torques) on the four arms and body. For the 
symmetric swimmer, the torque and force in the transversal direction vanish 
by symmetry. The swimming equation is the condition that the force in the 
body-direction vanishes. This force depends linearly on the known rate of 
change of the controls, 9 and the unknown the velocity X of the "body-rod" . 
It gives a linear equation for the velocity. 

For slender arms and body, the forces are given by Cox j5] slender body 

1 Multimedia simulations can be viewed in http : //physics . technion . ac . il/~avron 
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theory: The element of force, dF(s), acting on a segment of length ds located 
at the point s on the slender body is given by 

dF(s) — k(t(t ■ v) - 2v) ds, k = ^ (2.1) 

In k 

where t(s) is a unit tangent vector to the slender-body at s and v(s) its 
velocity there, /x is the viscosity and k the slenderness (the ratio of length 
to diameter). 

The force on the a-th arm depends linearly on the velocities of the controls 
6j and swimming velocity X. For example, the force component in the x- 
direction on the a-th arm takes the form 

FZ = fZ j O j + f: x X (2-2) 
where /„ ■ are functions of the controls, given by elementary integrals 

/** = A;(cos 2 9 a - 2) I ds , /* = 2k sin 6 j I s ds (2.3) 
Jo Jo 

Similar equations hold for the left arm and the body. The requirement that 
the total force on the swimmer vanishes gives a linear relation between the 
variation of the controls and the displacement dx 

dx 

dx= — = -a(£, 77) (d£i + d£ 2 ), & = cos 9 V (2.4) 

where 

As one expects, the body is just a "dead weight" and a trim swimmer with 
i] = is best. 



The curvature 

The notation dx stresses that the differential displacement does not integrate 
to a function of the controls, x(£). This is the essence of swimming: x fails 
to return to its original value with £. For this reason, swimming is best 
captured not by the differential one-form dx but by the differential two-form 
T = ddx 

T = dAx = 2a 2 (£, 77) (6 - 6) d& A d& (2.6) 
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Figure 3: The swimming curvature (measured in the dissipation metric) on 
the surface of controls for the symmetric Purcell swimmer. The closed curve 
shows the optimal stroke. A bona fide optimizer exists since the swimming 
curvature vanishes on the boundary 9 a = tt, . 



and a(£) is the rational function given in Eq. (2.5). T is commonly known 



as the curvature [H] and its surface integral for a region enclosed by a curve 
7 gives, by Stokes, the distance covered in one stroke. It gives complete 
information on both the direction of swimming and distance. 

The total curvature associated with the full square of shape space, |£j | < 1 
is by symmetry. The total positive curvature associated with the triangular 
half of the square, £i < £2 is 0.274. This means that swimmer can swim, at 
most, about a quarter of its arms length in a single stroke. 

T is a differential two form and as such is assigned a numerical value only 
when in comes with a region of integration. By itself, it has no numerical 
value. To say that the curvature is large at a given point in shape, or control, 
space requires fixing some a-priori measure. For example one can pick the 
flat measure for £ in which case a 2 (£,7/)(£i — £2) gives numerical values for 
the curvature. However, one can pick instead the flat measure for 9 one gets 
a different function. A natural measure on shape space is determined by the 
dissipation. We turn to it now. 
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The metric in shape space 

The power of swimming at low Reynolds numbers is quadratic in the driving: 
QjkOjOk where gjk(0) is a function on shape space and we use the summa- 
tion convention where repeated indices are summed over. This suggests the 
natural metric in shape space is, in either coordinate systems, 

g jk (6) dOjdOk = g jk (0 d£,-d& (2.7) 

In particular, the associated form is 



y/det g{6) d6 x A d0 2 = V det 9(0 d & A d & ( 2 - 8 ) 
The curvature can now be assigned a natural numerical value 



2a 2 (£,r,) ^jl= = 2a 2 (e,r / )^=%=sing 1 sing 2 (2.9) 
V det 9(0 V d etg{6) 

Each arm of the symmetric Purcell swimmer dissipate energy at the rate 

- dF(s)-vds — —k ((t ■ v) 2 - 2v ■ v) ds (2.10) 
Jo Jo 

= - k [ (X 2 - 2s 2 e] cos 2 9j - 2(s9j sin 6 j - X) 2 ) ds 
Jo 

= (:5J' 2 • -lit] ■ (iC r r) 
And the total energy dissipation by the arms is evidently 

— (Qx 2 + 2(61 + 01) + 6(6 + (2-11) 

In a body-less swimmer, 7] = 0, this is also the total dissipation, and we 
consider this case from now on. Since we are interested in the metric up to 
units, we shall henceforth setj^] 4k£ 3 = 3. Plugging the swimming equation, 
Eq. (2.4) gives g: 

/n\ t<-*\( 5-2£f + £ 2 sin Q x sin 6 2 \ t n 
9(0) = a(e, 0) ( sin di 2 in Q l 5 _ ^ + e ), = cos 9 3 (2.12) 



and a(6 0) is given in Eq. (2.5). In particular, g{6) is a smooth function on 
shape space while <?(£) is singular at the boundaries. One can now meaning- 
fully plot the curvature which is shown in Fig. [3] 



2 An alternate natural normalization is one that preserves the area Air 2 



7 



The optimal stroke 

Efficient swimming covers the largest distance for given energy resource and 
at a given speecQ Alternatively, it minimizes the energy needed for covering 
a given distance at a given speecQ Fixing the speed for a given distance is 
equivalent to fixing the time r. In this formulation the variational problem 
takes the form of a problem in Lagrangian mechanics of minimizing the action 

T 

J 9jk{0) djOkdt + q J dx (2.13) 



where q is a Lagrange multiplier, dx is given in Eq. (2.4). This can be 
interpreted as a motion of a a charged particle on a curved surface in an 
external magnetic field [2]. Conservation of energy then says that the solution 
has constant speed (in the metric g). For a closed path, the kinetic term is 
then proportional to the length of the path and the constraint is the flux 
enclosed by it. Thus the variational problem can be rephrase geometrically 
as the "isoperimetric problem" : Find the shortest path that encloses the 
most flux. 

The charged particle moves on a curved surface. How does this surface 
look like? From the dissipation metric we can calculate, using Brioschi for- 
mula [32], the gaussian curvature K (not to be confused with JF) of the 
surface. A plot of it is given in Fig. |4| 

Inspection of Fig. [3] suggest that pretty good strokes are those that enclose 
only one sign of the curvature T . The actual optimal stroke can only be found 
numerically. It is plotted in figure [3} The efficiency for this stroke is about 
the same as the efficiency of the Purcell's swimmer for rectangular strokes of 
[3], but less than the optimally efficient strokes found in [16J. 

Cox theory does not allow the arms to get too close. How close they are 
allowed to get depends on the slenderness k. The smallest angle allowed 59 
must be such that (259) log k ^> 1. As the optimal stroke gets quite close 
to the boundary, with 59 ~ 0.1 radian it can be taken seriously only for 
sufficiently slender bodies with log k ^> 5, which is huge. The optimal stroke 
is therefore more of mathematical than physical interest. One can use a refine 



3 Onc needs to constrain the average speed since one can always make the dissipation 
arbitrarily small by swimming more slowly. 

4 The distance is assumed to be large compared with a single stroke distance. The 
number of strokes is not given a-priori. 
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slender body approximation by taking high order terms in Cox's expansion 
for the force. This will leave the structure without changes, but will made 
Eq. |2.6 and Eq. 2.12 much more complicated. 

From a mathematical point of view it is actually quite remarkable that 
a minimizer exists. By this we mean that the optimal stroke does not hit 
the boundary of shape space |^| = 1 where Cox theory is squeezed out 
of existence. This can be seen from the following argument. Inspection 
of Eq. (2.9) and Eq. (2.12) shows that the curvature vanishes linearly near 
the boundary of shape space (this is most easily seen in the 9 coordinates). 
Suppose now that the optimal path ran along the boundary. Shifting the 
path a distance e away from the boundary would shorten it linearly in e 
while the change in the flux integral will be only quadratic. This shows that 
the path that hits the boundary can not be a minimizer. 



3 Purcell swimmer 

Purcell swimmer can move in either direction in the plane and can also rotate. 
Since the Euclidean group is not Abelian (rotations and translations do not 
commute) the notion of "swimming curvature" that proved to be so useful in 
the Abelian case needs to be modified. As we shall explain, landscape figures 
can be used to give qualitative geometric understanding of the swimming and 
in particular can be used to answer Purcell question "What will determine 
the direction this swimmer will swim?". However, unlike the Abelian case, 
the swimming curvature does not give full quantitative information on the 
swimming and one can not avoid solving a system of differential equations 
in this case if one is interested in quantitative details. 

The location and orientation of the swimmer (in the Lab frame) shall be 
denoted by the triplet x a where x° = <fr is the orientation of the swimmer, 
see Fig. [TJ and x 1 ' 2 are cartesian coordinates of the center of the "body" 
We use super-indices and Greek letters to designate the response while lower 
indices and Roman characters designate the controls \6j\ < n. 

Linear response versus a gauge theory 

The common approach to low Reynolds numbers swimming is to write the 
equations of motion in a fixed, lab frame. We first review this and then 

5 A11 distances are dimcnsionless being measured in units of arm length i. 
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describe an alternate approach where the equations of motions are written 
in a frame that instantaneously coincides with the swimmer. 

By general principles of low Reynolds numbers hydrodynamics there is a 
linear relations between the change in the controlaj d6 and the response dx 



dx a = AJd9 :j , (3.1) 

(summation over repeated indices implied.) Note that j = 1,2 since there 
are two controls, while a = 0, 1, 2 for the three responses. 

The response coefficients A are functions of both the control coordinates 
8k and the location coordinates x 13 of the swimmer in the Lab. However, 
in a homogeneous medium it is clear that A® can only be a function of the 
orientation x° = 0. Moreover, in an isotropic medium it can only dependent 
on the orientation through 



A?((j),9) = R al3 ((j))A?(6); R af3 (0) = cos0 -sin0 (3.2) 




In the Lab frame, the nature of the solution of the differential equations 
is obscured by the fact that one can not determine dx independently for 
different points on the stroke (because of the dependence on <p = x°). 

The coefficients A? may be viewed as the transport coefficients in a rest 
frame that instantaneously coincides with the swimmer. They play a key 
role in the geometric picture that we shall now describe. In the frame of the 
swimmer one has 

dy a = A«d6 j (3.3) 

which is an equation that is fully determined by the controls. The price one 
pays is that the dy coordinates cannot be simply added to calculate the total 
change in a stroke 7, since one has to consider the changes in the reference 
frame as well. In order to do that, dy must be viewed as (infinitesimal) 
elements of the Euclidean group 

cos y° sin y° 
E(y a ) = j -siny cosy y 2 j (3.4) 




6 In this section we use the convention the 9\ increases counterclockwise and #2 clock- 
wise. 
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The composition of dy along a stroke 7 is a matrix multiplication 

E(j) = l[E(d y <*(e)) (3.5) 

6>S7 

The product is, of course, non commutative. We denote generators of trans- 
lations and rotations by 






e° = -1 , e 1 = I I , e 2 = | 1 | (3.6) 
\ / 

They satisfy the Lie algebra 

[e°,e a ] = -e 0af> e f> } [e\e 2 ]=0 (3.7) 

where e a ^ is the completely anti-symmetric tensor. One can write Eq. (3.3) 
concisely matrix equation 

dy = AjdOj, dy = y a e a , A j = A«e a (3.8) 

where dy and Aj are 3x3 matrices (summation over repeated indices implied). 

The swimming curvatures 

Once the (six) transport coefficients A" are known, one can, in principle, 
simply integrate the system of three, first order, non-linear ordinary differ- 



ential equations, Eq. (3.1). This can normally be done only numerically. 
Numerical integration is practical and useful, but not directly insightful. We 
want to describe tools that allow for a qualitative understating swimming in 
the plane without actually solving any differential equation. 

Low Reynolds numbers swimmers perform lots of mutually cancelling ma- 
neuvers with a small net effect. The swimming curvature measure only what 
fails to cancel for infinitesimal strokes. Since reversing a loop reverses the 
response, it is natural to expect, that Sy a for a closed (square) loop is propor- 



tional to the area form. Integrating Eq. (3.8) around a closed infinitesimal 
loop gives 

5y = TdQ x A d6 2 (3.9) 

where 

T = dxA 2 — d 2 Ax - [Ai, A 2 ], T = T a e a , d 5 = A (3.10) 
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T and Sy are 3x3 matrices. T has the structure of curvature of a non-abelian 
gauge field [12J. In coordinates, this reads 



T a = d x A% - d 2 A1 + e 0a/3 (A\A P 2 - A°a() , (3.11) 
In the Lab coordinates one has, of course, 

Sx a = T^dQx A d02, T a = R afB (0)^(9), (3.12) 
The curvature is Abelian when the commutator vanishes. This is the case 



in Eq. (2.9) and it is also the case for the rotational curvature. The Abelian 
curvature gives full information on the swimming of finite stroke by simple 
application of Stokes formula. This is, unfortunately, not the case in the 
non- Abelian case. One can not reconstruct the translational motion of a 
large stroke from the infinitesimal closed strokes Sy because Stokes theorem 
only works for commutative coordinates and Sy are not. 

For Purcell swimmer, T although explicit, is rather complicated. Since 
the dissipation metric is complicated too, we give two plots of T: Figs. (|5j 
8p ) give the curvature relative to the flat measure on (61,62), and describe 



how far the swimmer swims for small strokes. In Figs. (|6j 10|11 ) the curva- 
ture is plotted relative to the dissipation measure and it displays the energy 
efficiency of strokes. 

The equations of motion 

Using Cox theory, in a manner analogous to what was done in for the sym- 
metric swimmer, one can calculate explicitly the force (and torque) on the 
a-th rod in the form 

K = m + ft^ (3-13) 
where are explicit and relatively simple functions of the controls (compare 



with Eq. (2.3)). The swimming equation are then given by 



£ F " = ( £ ft) °i + ( E " J ±p = ( 3 - 14 ) 

a \ a J \ a / 

This reduces the problem of finding the connections A to a problem in linear 
algebra. Formally 

a'=(yjA (yj:] (3.15) 
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where the bracket on the left is interpreted as a 3 x 3 matrix, with entries 
a,f3, and the inverse means an inverse in the sense of matrices. Although 
this is an inverse of only a 3 x 3 matrix the resulting expressions are not 
very insightful. We spare the reader this ugliness which is best done using a 
computer program. 



Symmetries 

Picking the center point of the body as the reference fiducial point is, in the 
terminology of Wilczek and Shapere [IB] a choice of gauge. This particular 
choice is nice because it implies symmetries of the connection A [THl [3]. 
Observe first that the interchange (61,62) — > (—62,-61) corresponds to a 



rotation of the swimmer by ir. Plugging this in Eq. (3.2) one finds 



+Al(-6 2 ,-6 1 ), (3 = 0; 
- A 2 (— 62, — 61) , otherwise . 

This relates the two half of the square divided by the diagonal 61 + 6 2 = 0. 

A second symmetry comes from the interchange (6\, 62) — > (62, 61) corre- 
sponding to the reflection of the swimmer around the central vertical of the 
middle link. Some reflection shows then that 

aP/q n \ _ f + A 2^2,6l), (3 = 2; 

M9l > 92) - [-4(62,6,), otherwise. (3 - 17) 

This relates the two halves of the square divided by the diagonal 6\ = 62- 

The symmetries can be combined to yield the result that A and A 2 are 
anti-symmetric and A 1 is symmetric under inversion 

A%6) = -A%-9), A)(6) = A)(-6), A)(6) = -A)(-6) (3.18) 



Rotations 

The rotational motion of Purcell swimmer, in any finite stroke, is fully cap- 
tured by the Abelian curvature 

F° = F° = d x A° 2 - d 2 A\ (3.19) 

This reflects the fact that rotations in the plane are commutative. 
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The symmetry of Eq. (3.16) implies 



(d 2 A° 1 )(e 1 ,e 2 ) = (d 1 A° 2 )(-e 2 ,-e 1 ) (3.20) 

and this says that J-°is ani- symmetric under reflection in the diagonal Q\ + 



9 2 = 0. Similarly, Eq. (3.21) implies 



(d 2 A° 1 )(e 1 ,e 2 ) = -(d 1 A° 2 )(e 2 ,e 1 ), (3.21) 

and this says that JF° is symmetric about the line Q\ = Q 2 . Fig. ^ is a plot 
of the curvature and it clearly has the requisite symmetries. 

The total curvature associated with the full square of shape space vanishes 
(by symmetry). For rj — 2, one can see in Fig. [5] three positive islands 
surrounded by three negative lakes. The total curvature associated with the 
three islands is quite small, about 0.1. This means that Purcell swimmer 
with rj = 2 turns only a small fraction of a circle in any full stroke. 



Translation 



The curvatures corresponding to the two translations of a swimmer with rj = 
0.75 are shown in Figs. ( 8|9|l0|ll )(here we use r] = 0.75 for comparison with 
[16]). The symmetries of the figures are a consequence of Eqs. ( 3.16|3.21 ). 
Form the first we have 



(d 2 Af)(6 u 8 2 ) = -(d 1 A%)(-9 2 ,-9 l ),P= 1,2 



(3.22) 



which implies that JF 1,2 are symmetric under reflection in the diagonal Q\ + 
9 2 = 0. Similarly, from the Eq. (3.21) we have 



(d 2 Af)(e 1 ,e 2 



-(d 1 AD(e 2 ,e l ), (3 = 2; 

-(d 1 Al)(9 2 ,9 1 ), (3 = 1. 



(3.23) 



This says that T 1 symmetric and T 2 anti- symmetric under reflection in the 
diagonal Q\ = 9 2 . 

The curvatures for the translations is non-Abelian and can not be used 
to calculate the swimming distances for finite strokes because the Stokes 
theorem fails. 
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4 Qualitative analysis of swimming 



When is a stroke small? 

The landscape figures for the translational curvatures provide precise infor- 
mation on the swimming distance for infinitesimal strokes. They are then 
also useful to characterize small strokes. The question is how small is small? 
For a stroke of size e, the controls are of size 58 = 0(e) and the swimming 
distance measures by the curvature is 0(jFe 2 ). The error in this has term^] 
of the form 0(AFe 3 ). This suggest that the relative error in the swimming 
distance as measured a finite stroke is of the order O(Ae). Hence, a stroke 
is small provided \As\ <C 1. Clearly, a Purcell swimmer swims substantially 
less than an arm length as the arm moves. This says that \A\ <C 1 and so 
strokes of the order of a radian can be viewed as small strokes. 

A radian is the scale of the structures in the landscape of the figures of 
the curvature. This means that the landscape carries qualitative information 
about the swimming of moderate strokes. 

x versus y 

The x-curvature is symmetric under inversion 



Since both A and A 2 are antisymmetric under inversion, one sees that the 
non-Abelian part of the x-curvature is of order 0(6 2 ) near the origin. The x- 
translational curvature, which is non-zero near the origin, is almost Abelian 
for small strokes. 

The y-curvature, in contrast, is anti-symmetric under inversion 



and so vanishes linearly at the origin. The non-Abelian part is also anti- 
symmetric under inversion, and it too vanishes linearly. The y-curvature is 
therefore not approximately Abelian for small strokes, but it is small. 

7 There are also terms of the form the form 0(A 3 e 3 ). 




(4.1) 




(4.2) 
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Which way does a swimmer swim? 

The swimming direction can be easily determined for those strokes that live 
in a region where the translational curvature has a fixed sign. This answers 
Purcell's question for many strokes. Strokes the enclose both signs of the 
curvature are subtle. 



Subtle swimmers 



Purcell's swimmer can reverse its direction of propagation by increasing the 
stroke amplitude [3]. This can be seen from the landscape diagram, Fig. (j8|: 
Small square strokes near the origin sample only slightly negative curvature. 
As the stroke amplitude increases the square gets larger and begins to sam- 
ple regions where the curvature has the opposite sign, eventually sampling 
regions with substantial positive curvature. 



Optimal distance strokes 

The curvature landscapes are useful when one wants to search for optimal 
strokes as they provide an initial guess for the stroke. (This initial guess can 
then be improved by standard optimization numerical methods.) 

For example, Tarn and Hosi [TH] looked for strokes that cover the largest 
possible distance. For strokes near the origin, a local optimizer is the stroke 
that bounds the approximate square blue region in Fig. [8] (in this case r] = 
0.75). 



Efficient strokes 



The curvature normalized by dissipation, Fig. 10 gives a guide for finding 
efficient small strokes. Caution must be made, since while the displacement 
can be approximated from the surface area, the energy dissipation is pro- 
portional to the stroke's length and not the stroke's area. In regimes where 



the Gaussian curvature of the dissipation (Fig. 12) is positive - it is possible 
to have strokes with small length which bounds large area. In the case of 
Purcell's swimmer, this suggest two possible regimes: around the origin and 



the positive curvature island in the upper left (lower right) corner of Fig. 12 
The optimizer near the origin is the optimal stroke found in [16], while the 
optimizer in the upper left corner [1 



to the values of -, = in Fig. 

y/detg(0) 



ID 



- although more efficient (pay attention 
is of mathematical interest only, since 
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it is near the boundary, where the first order slender body approximation 
Eq.( |2.1 ) is relevant only for extremely slender bodies. 
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Figure 4: The Gaussian curvature on the surface of controls induced by the 
dissipation metric. 
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Figure 5: The curvature for rotation for Purcell's three linked swimmer with 
rj = 2, plotted with the flat measure on 2 ) 
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Figure 6: The rotation curvature for Purcell's three linked swimmer with 
7] = 2 plotted using the measure induced by dissipation 




Figure 7: The failure of Stokes theorem in the non-commutative case: The 
integration on the adjacent segments traversed in opposite senses do not 
cancel in the non-commutative case. 
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Figure 8: The landscape of the x-curvature for Purcell's swimmer with 77 = 
0.75 shown with the the Tarn and Hosoi optimal distance stroke [TB] . The 
curvature is given relative to the flat measure in (9±, 9 2 ). 
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Figure 9: The landscape of the y-curvature for Purcell's swimmer with 77 = 
0.75. The curvature is given relative to the flat measure in (#1, 9 2). 
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Figure 11: The y-curvature relative to the dissipation metric. Pay attention 
to the large values in the scale. 
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